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In a recent paper by Link, it was pointed out that the standard picture of the neutron star core 
composed of a mixture of a neutron superfluid and a proton type-II superconductor is inconsistent 
with observations of a long period precession in isolated pulsars. In the following we will show that 
an appropriate treatment of the interacting two-component superfluid (made of neutron and proton 
Cooper pairs), when the structure of proton vortices is strongly modifled, may dramatically change 
the standard picture, resulting in a type-I superconductor. In this case the magnetic fleld is expelled 
from the superconducting regions of the neutron star leading to the formation of the intermediate 
state when alternating domains of superconducting matter and normal matter coexist. 



PACS numbers: 

The conventional picture of a neutron star is that the 
extremely dense interior is mainly composed of neutrons, 
with a small amount of protons and electrons in beta 
equilibrium. The neutrons form '^P2 Cooper pairs and 
Bose condense to a superfluid state, while the protons 
form Cooper pairs and Bose condense, as well, to give 
a superconductor (see e.g. for a review). It is generally 
believed that the proton superfluid is a type-II supercon- 
ductor, which means that it supports a stable lattice of 
magnetic flux tubes in the presence of a magnetic field. In 
addition, the rotation of a neutron star causes a lattice 
of quantized vortices to form in the superfluid neutron 
state, similar to the observed vortices that form when 
superfluid He is rotated fast enough. In a recent paper 
by Link , it was pointed out that the precession of the 
neutron star hints that this picture may not necessarily 
be correct. In particular, Link states that the observed 
precession of a neutron star does not allow the proton 
magnetic flux tubes and neutron vortex lattice to exist 
simultaneously, due to the fact that the axis of rotation 
and the axis of the magnetic field are not aligned and 
the fact that these two different vortices interact quite 
strongly. Furthermore, Link suggests that the conven- 
tional picture of a neutron star as a type-II supercon- 
ductor may have to be reconsidered. One should remark 
here that the conventional picture of type-II supercon- 
ductivity follows from the standard analysis when only a 
single proton field is considered. As we shall demonstrate 
in this letter, if one takes into account that the Cooper 
pairs of neutrons are also present in the system and that 
they interact strongly with the proton Cooper pairs, the 
superconductor may in fact be type-I and exhibit the 
Meissner effect (total expulsion of an external magnetic 
field) , contrary to the picture that is obtained when only 
the proton Cooper pair condensate is accounted for. This 
would support the suggestion made by Link Q that neu- 
tron stars may in fact be type-I superconductors with the 
superconducting region not carrying any magnetic flux. 

The core of a neutron star is a mixture of neutron and 



proton superfluids, as discussed above. In the presence 
of a magnetic field, it is well known that the type-II pro- 
ton superfluid may form magnetic flux tubes. Inside the 
core of these vortices, the proton condensate vanishes, 
and the core is filled with normal protons resulting in 
the restoration of the broken U{1)em symmetry. If the 
accepted estimates of the proton correlation length and 
the London penetration depth are used, then the distant 
proton vortices repel each other leading to formation of a 
stable vortex lattice. This is the standard picture realized 
in conventional type-II superconductors. However, there 
are many situations where this picture will be qualita- 
tively modified. For example, if a second field or compo- 
nent is added, such that there is an approximate SU (2) 
symmetry between the original and the second fields, it 
may be energetically favorable for the second field to con- 
dense inside the vortex core pj, resulting in a different 
pattern of the vortex-vortex interaction. This behavior 
is known to occur in various systems: cosmic strings, 
high Tc superconductors, Bose-Einstein condensates, su- 
perfluid ^He, and high baryon density quark matter (see 
Refs. iSiiH). 

In the case considered in this letter, we have a situa- 
tion where there are two condensates, proton and neutron 
Cooper pairs, both of which are nonzero in the bulk of 
the matter. In what follows we shall argue that if the in- 
teractions between the proton and neutron Cooper pairs 
at small momentum are approximately equal (a precise 
condition of "approximately" will be derived below) , the 
vortex-vortex interaction will be modified and the sys- 
tem will be a type-I superconductor with the magnetic 
field completely expelled from the superconducting re- 
gions. We believe that the approximate symmetry of pro- 
ton/neutron Cooper pair interactions at large distances 
is somewhat justified by the original isospin symmetry of 
bare protons and neutrons, however this symmetry is not 
exactly equivalent to the conventional isotopical SU{2) 
symmetry. If we consider a proton vortex (magnetic flux 
tube) in this case, the vortex structure is non-trivial, as 
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we will see below. The core of the proton vortex, where 
the proton superfluid density goes to zero, has a neutron 
superfluid density that is larger than at spatial infinity, 
far from the core. Moreover, the size of the vortex core 
and the asymptotic behavior of the proton condensate far 
from the core are also modified due to the additional neu- 
tron condensate that is present. The most important re- 
sult of these effects is that the interaction between distant 
proton vortices may be attractive in a physical region of 
parameter space leading to type-I behavior: destruction 
of the proton vortex lattice and expulsion of the magnetic 
flux from the superconducting region of the neutron star. 
We will now elaborate on the ideas outlined above. 

We start by considering the following effective Landau- 
Ginsburg free energy that describes a two component 
Bose condensed system. In our system, we have a pro- 
ton condensate described by the field ipi and a neutron 
condensate described by the field tp2- The tpi field with 
electric charge q (which is twice the fundamental pro- 
ton charge, q = 2|e|) interacts with the gauge field A, 
with B = V X A. The two dimensional free energy reads 
(we neglect the dependence on third direction along the 
vortex) : 

J Ziric he 

+ ^ + ni^iiMV'2p)], (1) 

where TOc = 2m and m is the mass of the nucleon. Here 
we have moved the effective mass difference of the proton 
and neutron Cooper pairs onto the interaction potential 
V . In the free energy given above, we have ignored the 
term coupling the proton and neutron superfluid veloc- 
ities, which gives rise to the Andreev-Bashkin effect 8], 
as it is not important in our discussions. Indeed, the 
relevant term in the free energy can be represented as 
~ / (Pxvi ■ V2, where vi and V2 are velocities of the su- 
perfluid components. For neutron stars which do not ro- 
tate (the case which is considered in this paper) V2 — 0, 
and the effect obviously vanishes. We expect that due 
to the small density of the neutron vortices (compared 
to the density of the proton vortices) the effect is still 
negligible for most of the flux tubes in a rotating star as 
well. The effect could be important only for a few of the 
flux tubes situated close to a neutron vortex core, where 
V2 strongly deviates from the constant value at interflux 
distance scales. 

We have also ignored the fact that the neutron con- 
densate has a non-trivial ^P2 order parameter as only the 
magnitude of the neutron condensate is relevant to the 
effect described below. The free energy Q only describes 
large distances and it does not describe the gap structure 
on the Fermi surfaces, only the superfluid component of 
the protons and neutrons. 

The free energy ^ is invariant under a C/(l)i x U{1)2 
symmetry associated with respective phase rotations of 



fields '01 and ip2, which corresponds to the conservation of 
the number of Cooper pairs for each species of particles. 
Moreover, we know that the free energy describes par- 
ticles interacting via the strong nuclear force and, there- 
fore, must be approximately invariant with respect to 
the SU{2) isospin symmetry. Therefore, we expect that 
while the Fermi surfaces for protons and neutrons are 
very different and the gap equations are very different, 
the interaction between different Cooper pairs at small 
momentum must not be very different (the asymmetry 
must be proportional to (m^ — m„)). This asymmetry 
is expressed in terms of different scattering lengths of 
Cooper pairs for each species. Thus, we assume, the in- 
teraction potential V can be approximately written as 
V'dV'iP, |V'2p) ~ C/(|^ip-)-|V'2p)- In reality this symme- 
try is explicitly slightly broken, and the potential V has 
a minimum at IV'iP — "ii |V'2p = '^2 where ni and n2 
are the proton and neutron Cooper pair densities. Hence 
in the ground state, \ < ipi > — rii, i — 1, 2, and both 
C/(I) symmetries are spontaneously broken. An impor- 
tant quantity for the analysis that follows will be the ratio 
of proton to neutron Cooper pair density, 7 = ni/n2. A 
typical value of 7 in the core of a neutron star is 5— 15%, 
thus, in our numerical estimates below we will often use 
the limit 7^1, though our qualitative results do not 
depend on this parameter |l4| . 

Now let's investigate the structure of proton vortices, 
which exist due to the spontaneous breaking of the U{l)i 
symmetry. Such vortices are characterized by the phase 
of the tpi field varying by an integer multiple of 27r as 
one traverses a contour around the core of the vortex. 
By continuity, the field ipi must vanish in the center of 
the vortex core. Up to this point, it has been assumed 
that the neutron order parameter ■02 will remain at its 
expectation value in the vicinity of the proton vortex. As 
we have already remarked, this is not the case in many 
similar systems, see e.g. |3, Ij, lu III ) • 

So, anticipating a non-trivial behavior of the neutron 
field 02, let's adopt the following cylindrically symmetric 
ansatz for the fields describing a proton vortex with a 
unit winding number: 

01 = V^/(r)e^^ V2 = V^ff(r), A = -^9 (2) 

q r 

where (r, 0) are the standard polar coordinates. Here we 
assume that the proton vortex is sufficiently far from any 
rotational neutron vortices, so that any variation of 02 
is solely due to the proton vortex. The functions /, g, 
and a obey the following boundary conditions: /(O) — 0, 
/(oo) = 1, 5'(0) = 0, g{oo) = 1, a(0) = 0, and a(oo) = 1. 
We see that the fields 0i and ^02 approach their expec- 
tation values at r = 00. 

We wish to find the asymptotic behavior of fields 0i , 
"02, and A far from the proton vortex core, as this will 
determine whether distant vortices repel or attract each 
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other. The asymptotic behavior can be found analyti- 
cally by expanding the fields defined in jSJ: 

f{r) = 1 + F(r), g{r) = 1 + G(r), a(r) = 1 - r5(r) (3) 

so that far away from the vortex core, F, G,rS <^ 1 and 
F,G,S^Oasr—^oo. This allows us to linearize far 
from the vortex core the equations of motion correspond- 
ing to the free energy to obtain: 
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(4) 
(5) 



S" + -5' 
r 

where the London penetration depth A ~ 
\J rricC^ / ^nq^ni. Here all derivatives are with re- 
spect to r, and the matrix M mixing the fields F and G 
is, 
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where the derivatives Vij = d'^V/{d\tpi\'^d\'ipj\-^) are eval- 
uated at IV'iP = rii. Here we assume that <^ F,G, 
i.e. the superconductor is not in the strong type-II 
regime (this is justified since we are only attempting 
to find the boundary between type-I and type-H super- 
conductivity). The solution to Eq. Q is known to be: 
S = ^Ki{r/X) where Ki is the modified Bessel func- 
tion and Ga is an arbitrary constant. The remaining 
equation Q can be solved by diagonalizing the mixing 
matrix M. In previous works the influence of the neutron 
condensate on the proton vortex was neglected, which 
formally amounts to setting the off-diagonal term V12 in 
M to 0. In that case, one can assume that the neutron 
field remains at its expectation value, i.e. G = 0, to ob- 
tain, F = GFK(){\/2r /£) where ^ — ^ j2mcn^Vi-Y is 
the correlation length of the proton superconductor and 
Kq is the modified Bessel function. It is estimated that 
A 80 fm and ^ ^ 30 fm which leads to k = Xji^ ^ Z 
for the Landau- Ginzburg parameter. As is known from 
conventional superconductors, if k > 1/a/2, distant vor- 
tices repel each other leading to type-II behavior. This 
is the standard picture of the proton superconductor in 
neutron stars that is widely accepted in the astrophysics 
community. 

However, the standard procedure described above is in- 
herently flawed since the system exhibits an approximate 
C/(2) symmetry, which forces approximate equality of sec- 
ond partial derivatives, V\\ ~ V22 ~ V12. This makes the 
mixing matrix M nearly degenerate. The general solu- 
tion to Eq. Q is: 



j = ^ GiKo{y/T^r)vi 



(7) 



where Vi and Vj are the eigenvalues and eigenvectors 
of matrix M, and Gi are constants to be calculated 



by matching to the solution of the original non-linear 
equations of motion. In the limit 7 = ni/n2 <C 1 and 
e = (V11V22 — ^iD/^i^ 1 0116 can estimate : 
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— , VI ~ (-1,7), V2~(l,l). (8) 



The physical meaning of solution Q is simple: there are 
two modes in our two component system. The first mode 
describes fluctuations of relative density (concentration) 
of two components and the second mode describes fluctu- 
ations of overall density of two components. Notice that 
1^1 V2, and hence the overall density mode has a much 
smaller correlation length than the concentration mode. 
Therefore, far from the vortex core, the contribution of 
the overall density mode can be neglected, and one can 
write: 
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The most important result of the above discussion is that 
the distance scale over which the proton and neutron con- 
densates tend to their expectation values near a proton 
vortex is of order ^/-^e - the correlation length of the con- 
centration mode. Since e <C 1, this distance scale can be 
much larger than the proton correlation length ^, which 
is typically assumed to be the radius of the proton vortex 
core. 

We have also verified the above results numerically 
by solving the equations of motions corresponding to 
with a particular choice of the approximately U{2) sym- 
metric interaction potential V. Our numerical results 
support the analytical calculations given above. Namely, 
we find that the magnitude of the neutron condensate 
is slightly increased in the vortex core, the radius of the 
magnetic flux tube is of order A and the radius of the 
proton vortex core is of order ^/\/e. 

Now that we know the approximate solution for the 
proton vortex, we will proceed to look at the interaction 
between two proton vortices that are widely separated. If 
the interaction between two vortices is repulsive, it is en- 
ergetically favorable for the superconductor to organize 
an Abrikosov vortex lattice with each vortex carrying a 
single magnetic flux quantum. As the magnetic field is 
increased, more vortices will appear in the material. This 
is classic type-II behavior. If the interaction between two 
vortices is attractive, it is energetically favorable for n 
vortices to coalesce and form a vortex of winding num- 
ber n, which is expelled from the sample. This is type- 
I behavior. Typically, the Landau- Ginzburg parameter 
K — X/£^ is introduced. In a conventional superconduc- 
tor, if K < l/\/2 then the superconductor is type-I and 
vortices attract. If k > 1/a/2 then vortices repel each 
other and the superconductor is type-II. As mentioned 
above, the typical value for a neutron star is k ~ 3, so 
we would naively expect that the proton superfluid is a 
type-II superconductor. 
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Now we will present three different calculations sup- 
porting our claim that for the typical parameters of a 
neutron star the proton superconductor may be type-I 
rather than type-II. First of all, we follow the method 
suggested originally in |0| to calculate the force between 
two widely separated vortices. The idea of this method 
is to model distant vortices as point sources in a free the- 
ory, which accurately describes the behavior of fields far 
from the vortex cores. The methods of were subse- 
quently applied in to the case similar to ours, the 
interaction of two widely separated vortices that have 
nontrivial core structure. Using the asymptotic field so- 
lutions found above, we follow the procedure of to 
obtain the following expression for the interaction energy 
per unit vortex length of two distant parallel vortices: 

U(d) c M!!^(c2x„(d/A) - CfKoiV2'ed/0) (10) 

where d — > oo is the separation between the two vortices. 
We see that if the first term in U dominates as d ^ oo 
then the potential is repulsive, otherwise, if the second 
term dominates the potential is attractive. We introduce 
the new dimensionless parameter Knp = VeA/^ into our 
description. In terms of this parameter, if Knp < 1/V2, 
then vortices attract each other and the superconduc- 
tor is type-I; otherwise, vortices repel each other and 
the superconductor is type-II. Therefore, our parameter 
n,ip = A/(5 = \/eA/^ should be considered as an effec- 
tive Landau- Ginzburg parameter, which determines the 
boundary between the type-I and type-II proton super- 
conductivity. Due to the importance and far reaching 
consequences of this result, we have also calculated the 
vortex- vortex interaction energy in a more direct way fol- 
lowing ^3 ; this calculation P| produced the same result 
(|10|l as the above procedure, therefore confirming our pic- 
ture. Our third check of the main result that for relatively 
small e the superconductor in the neutron stars may be, 
in fact, type-I is based on the macroscopical calculation 
of the critical magnetic fields. Usually one calculates the 
critical magnetic fields He and Hc2- We have calculated 
the critical magnetic fields He and Hc2 corresponding 
to the free energy ^ with a particular choice of V as 
a quadratic polynomial in IV'iP- The boundary between 
type-I and type-II superconductivity obtained using this 
procedure matches the results of our inter-vortex force 
calculation presented above. 

The most important consequence of this letter is that 
whether proton superconductor is type-I or II depends 
strongly on the magnitude of the SU (2) asymmetry pa- 
rameter e. Specifically, we find that the superconduc- 
tor is type-I when Knp = VeA/^ < l/\/2, and type-II 
otherwise. This result is quite generic, and not very 
sensitive to the specific details of the interaction poten- 
tial V. In particular, when e — > the superconduc- 
tor is type-I. The parameter e is not known precisely; 
the corresponding microscopical calculation would re- 



quire the analysis of the scattering lengths (amplitudes 
at small momentum) of Cooper pairs for different species. 
We can roughly estimate this parameter as being re- 
lated to the original SU (2) isospin symmetry breaking 
e ^ (m^ — mu)I^QCD ~ 10~^. If this is assumed to be 
the value of e, we estimate K„p = V^^/C 0.3 < l/\/2, 
which corresponds to a type-I superconductor. From 
these crude estimates, we see that it is very likely that 
neutron stars are type-I superconductors with the super- 
conducting region devoid of any magnetic flux, as was 
originally suggested in to resolve the inconsistency 
with observations of long period precession in isolated 
pulsars. If this is the case, some of the explanations of 
glitches [ISi (sudden changes of the neutron star's rota- 
tional frequency, see ^3 fo^' the original explanation of 
glitches), which assumes a type-II proton superconduc- 
tor, have to be reconsidered. Type-I superconductivity 
does not imply total expulsion of the magnetic field: the 
core structure could be composed of alternating domains 
of superconducting matter and normal matter 
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